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Abstract 

The algebraic structures underlying quasi-exact solvability for spin 
1/2 Hamiltonians in one dimension are studied in detail. Necessary 
and sufficient conditions for a matrix second-order differential opera- 
tor preserving a space of wave functions with polynomial components 
to be equivalent to a Schrodinger operator are found. Systematic sim- 
plifications of these conditions are analyzed, and are then applied to 
the construction of several new examples of multi-parameter QES spin 
1/2 Hamiltonians in one dimension. 
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1 Introduction 



Symmetries have traditionally played an essential role in quantum mechan- 
ics. For a few remarkable Hamiltonians, the knowledge of enough symmetries 
leads to a complete characterization of the spectrum by algebraic methods, 



14] . In general, however, the spectrum of an arbitrary Hamiltonian cannot 
be calculated analytically. During the last decade, a remarkable intermedi- 
ate class of quasi- exactly solvable (QES) spectral problems was introduced, 
for which a finite part of the spectrum can be computed by purely algebraic 



methods, |18|1 , [pljl , JIB] . The key feature in the latter class of spectral prob- 
lems is that the Hamiltonian H is expressible as a quadratic combination 
of the generators of a finite-dimensional Lie algebra q of first order differen- 
tial operators preserving a finite-dimensional module of smooth functions N '. 
Thus, H restricts to a linear transformation in the finite-dimensional vec- 
tor space A/", and therefore part of its spectrum can be computed by matrix 
eigenvalue methods. Appropriate boundary conditions must be imposed so 
that the eigenfunctions thus obtained qualify as physical wave functions, as 
e.g. square integrability if they represent bound states of the system, ||. 

These ideas, originally introduced for scalar Hamiltonians describing spin- 
less particles, can be generalized to include particles with spin. The first 



step in this direction was taken by Shifman and Turbiner, |jT? |, using the 



fact that a Hamiltonian for a spin 1/2 particle in d spatial dimensions can 
be constructed from a Lie superalgebra of first order differential operators 
in d ordinary (commuting) variables and one Grassmann (anticommuting) 
variable. Alternatively, JTJ, 2 x 2 matrices (or N x N matrices for parti- 
cles of arbitrary spin, ||) can be used to represent the Grassmann variable. 
However, in stark contrast with the scalar case, very few examples of ma- 
trix QES Schrodinger operators have been found thus far, |jl7l . There are 
two important conceptual reasons for this fact. First, the algebraic struc- 
tures underlying partial integrability in the matrix case are richer and less 
understood than in the scalar case. For one thing, as mentioned before, for 
matrix Hamiltonians Lie superalgebras of matrix differential operators nat- 
urally come into play, whereas in the scalar case only Lie algebras need be 
considered. Moreover, as we shall explain in Section |3], one even has to go 
beyond Lie superalgebras of matrix differential operators in order to explain 
quasi-exact solvability in the matrix case, [0. Secondly, 0, ||, every 
scalar second order differential operator in one dimension can be transformed 
into a Schrodinger operator of the form —d 2 x + V(x) by a suitable change of 
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the independent variable x and a local rescaling of the wave function. For 
matrix differential operators, the analogue of this result — V(x) being now a 
Hermitian matrix of smooth functions — is no longer true unless the operator 
satisfies quite stringent conditions, as we shall see in detail in Section |J 

The aim of this paper is to achieve a better theoretical understanding 
of quasi-exact solvability in the matrix case, which will enable us to con- 
struct new examples of matrix QES Schrodinger operators. To this end, in 
Sections |] and |3] we study the algebraic properties of certain algebras of ma- 
trix QES operators, reviewing the literature on the subject and obtaining 
several new results as well. In particular, we give a self-contained proof of 
the characterization of the class of QES matrix differential operators preserv- 
ing a finite-dimensional space of wave functions with polynomial components 
stated by Turbiner, [p0| , and Brihaye et al, [Q, 0. For the important partic- 
ular case of spin 1/2 particles, we derive in Section ^ necessary and sufficient 
conditions for a QES operator to be equivalent to a non-trivial Schrodinger 
operator. These conditions turn out to be too complicated to be solved in 
full generality, and so in Sections [5] and |6| we introduce some key simplifi- 
cations that will prove very useful in the task of finding explicit examples. 
Finally, the previous results are applied in Section ^ to the construction of 
several new examples of multi-parameter QES spin 1/2 Hamiltonians in one 
dimension. 



2 Scalar QES Operators 

We start with the scalar case, introducing the basic concepts and definitions 
and stating two theorems for the one-dimensional case which will play an 
important role in what follows. Since the results of this section are fairly 
standard, we will skip many details and all the proofs, referring the reader 
to the review articles [[RJ and [[/]] for an in-depth study. 

Let M denote an open subset of lR d , and let 2) 1 (M) be the Lie algebra of 
first order differential operators 

x = E ?V)-!n + iW> z = • • • ' *") e M ' 

acting on C°°(M), the Lie bracket being defined as the usual commutator 
between operators: 

[X,Y] = XY-YX, X,YeV\M). 
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Definition 2.1 A finite- dimensional Lie subalgebra g of D (M) is called 
quasi-exactly solvable (QES) if it preserves a finite- dimensional module N C 
C°°(M). A differential operator T is QES if it lies in the universal enveloping 
algebra 11(g) of a QES Lie algebra g. 

In general, quasi-exact solvability of a given differential operator T cannot 
be ascertained a priori. Therefore, the procedure usually followed consists 
in classifying QES Lie algebras modulo a suitable equivalence relation, and 
then using the canonical forms in the classification thus obtained to construct 
QES operators. 

Definition 2.2 Two differential operators T(z) and T(z) are equivalent if 

they are related by a change of the independent variables 



z = <p(z) (2.1) 
and a local scale transformation by a non-vanishing function U(z), i.e. 

T(z) = U(z)T(z)U- 1 (z). (2.2) 



The corresponding notion of equivalence for QES algebras follows di- 
rectly, i.e. two QES Lie algebras g and g are equivalent if their elements 
can be mapped into each other by a fixed transformation ( [2.1|) -( f2~2|) . Their 
associated finite-dimensional modules M and A/* are then related by 

j\f = U -Af, (2.3) 

the functions being expressed in the appropriate coordinates. The local clas- 
sification of finite-dimensional QES Lie algebras under the above notion of 
equivalence has already been completed for the case of one and two (real or 
complex) variables. Here we shall need only the one- dimensional case, fl3| . 

Theorem 2.3 Every (non-singular) QES Lie algebra in one (real or com- 
plex) variable is locally equivalent to a subalgebra of one of the Lie algebras 

Q n = Span{<9 2 , zd z , z 2 d z - nz, 1}, (2.4) 

where n G N. The associated g n -module is Af n = V n , the space of polynomials 
of degree at most n. 
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(The two-dimensional case, which is considerably more complicated but is 
not needed for the sequel, is discussed in 0, 0, and ||.) 

According to the previous theorem, every one-dimensional (scalar) QES 
differential operator T is locally equivalent to an operator T G U(Q n ) pre- 
serving V n for a suitable n. A partial converse of the latter result follows 
from the following remarkable theorem due to Turbiner, |R5|: 



Theorem 2.4 LetT^ be a k-th order linear differential operator preserving 
V n . We then have: 

i) Ifn>k, then may be represented by a k-th degree polynomial in 
the operators 

J+ = z 2 d z -nz } J° = z d z --, J- = d Z} (2.5) 

ii) If k > n, then = T <9™ +1 + T, where T is a linear differential 
operator of order k — n — 1, and T is a linear differential operator of 
order at most n satisfying i). 

The operators { J+, J°, J~} defined above span a QES Lie algebra Q n isomor- 
phic to s[(2), and the Lie algebras g n in Theorem [O] are simply a central 
extension by the constant functions of the corresponding g n . 



3 Algebraic Properties of PVSP Operators 

In the last section we have seen that every scalar QES scalar differential 
operator in one variable is essentially (up to equivalence) a polynomial in 
the generators of a Lie algebra g n preserving V n (for suitable n). When 
working with vector-valued wave functions, the natural generalization of V n 
is the polynomial vector space V nit ... >nN = V ni © • • ■ © V nN1 with elements 
ty(z) = fyi(z), . . . ,^N(z)y such that each component ipi is a polynomial of 
degree at most with complex coefficients. 

Definition 3.1 A N x N matrix differential operator T is called polynomial 
vector space preserving (PVSP) if it preserves V nu ... >nN = V ni © • • • © V nN 
for some non-negative integers n iy i = 1, . . . , N. 
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We will denote by the complex vector space of linear PVSP 

operators of order at most k preserving V nu ..., nN . Following U and 0, we 
will restrict ourselves in this paper to studying matrix PVSP differential 
operators. As we will be mainly concerned with spin 1/2 particles, the case 
N = 2 deserves special attention. 

3.1 Case N = 2 

Let n > A be non-negative integers, and consider the following set of matrix 
differential operators 



Q Q = z a (j~, Q Q = g Q (n, A)cr + , a = 0, . . . , A, 



with 



? a (n, A) = J] {zdz ~ n + A - k) (T z , (3.2) 

fc=i 

where we have adopted the convention that a product with its lower limit 
greater than the upper one is automatically 1, and a + = (er~)* = K ^ . 

It can be easily checked that the 6 + 2A operators in (3.1) (and also any 
polynomial thereof) preserve P„_A,n- We now introduce a Z 2 -grading in the 
set of 2 x 2 matrix differential operators D 2 x2 as follows: an operator 

T=( a b ) , a,b,c,de 2), 



c d. 



is said to be even if 6 = c = 0, and odd if a = d = 0. Therefore, the T's 
and J are even and the Q's and Q's odd. This grading, combined with the 
usual product (composition) of operators, endows T) 2X 2 with an associative 
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superalgebra structure. We can also construct a Lie superalgebra in D 2X 2 by 
defining a generalized Lie product by 

[A, B] s = AB — (-l) dcsAdcsB BA. (3.3) 

However, this product does not close within the vector space spanned by our 
operators (3.1), except for A = 0, 1. The explicit commutation relations are 
as follows, 1, [|: 



[J,T £ ] = 0, [J,Q a } = -^Q a , [J,Q a ] = ^Q a , 

[Qa, T £ ] = ( - a + | (1 + e))Q Q+e , [Q Q , T e ] = (a - | (1 - e))& 



{Qa, = {Qa, Qp} = 0, (3.4) 

where e = +, 0, — , and M a p is given, for a > (3, by 

A-a-l /3-1 

^/3= II {T° + J c -j-pP2)l[{T + J-k-{A-a)P 1 ) 

j=0 k=0 



with J c = A - 1 - J, and Pi = 1 - P 2 



'1 N 
0, 



As shown in [ffj, M aj g can be expressed in terms of T°, J, the identity, 
and the Casimir (for the even subalgebra) 



C = - l -(T+T- + T-T+) + T°T° = i f m(m + 2) , „. 
2 V y 4 1 n(n + 2) 



(3.5) 



where m = n — A, independently from n and the projectors Pi and P 2 . It 
can be readily verified that {Q a , Qp} gives a A-th order even differential 
operator, so the vector space spanned by the operators in (3.1) is not closed 
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under the Lie product ( |3.3| ) whenever A > 2. Moreover, it is not difficult 
to show that the Lie superalgebra Sa generated by the operators (3.1) is in 
this case infinite dimensional. Indeed, if we commute {Qa,Qo} = (^~) A 
with {Qo,Qa} = (T + ) A iteratively we obtain monomials in T + , T°, T~ of 
increasingly higher order. For A = 1 the underlying algebraic structure is 
the classical simple Lie superalgebra osp(2, 2), flSfl , [p0| , whereas for A = 



it is P)i ©sl(2), where b 1 is the 3-dimensional Heisenberg superalgebra. As 
remarked in |l| , in this latter case we can leave the grading aside and replace 
J = 1 by J = 03, ending up with the Lie algebra s[(2) © sl(2). 

Our next objective is to prove the analogous of Theorem |2]3] for PVSP 
operators preserving V n -A, n (first mentioned without proof by Turbiner for 
A = 1 in P0|] , an d subsequently by Brihaye and Kosinski for arbitrary A, 
Jl]]). It turns out that the operators (3.1) play the same role for the matrix 
case as the J's in ( |2.5p do for the scalar one. The proof relies on the next 
two Lemmas. 

Lemma 3.2 Let : V n — > V n -A be a k-th order linear differential opera- 
tor, with n > A. We then have: 

i) Ifn>k>A, then may be represented as 

T^ = J2q a (n,A)p a k _ A (Jn), (3-6) 

a=0 

where Pk-Ai^n) are polynomials of degree not higher than k — A in the 
operators J+, J° and J~ . 

u) J/A > jfe, then T&) = 0. 

Hi) If k > n, then T^ k ' = T d™ +1 + T, where T is a linear differential 
operator of order k — n — 1, and T is a linear differential operator of 
order at most n satisfying i) or ii). 

Proof, i). The essential point of the argument is that, as P„_a C V n , T^ k > 
preserves also V n , and so it must be a k-th degree polynomial in with 



suitable additional restrictions on its coefficients, according to Theorem \2A 
We will prove that 



TM=q (n,A)p° k _ A (J+,J° n ) + £ q a (n, A) p£_ A (J n °) 

Q=l 

+ q A (n, A) J' p k -A-i(J° n , J'), (3.7) 
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which clearly implies 

We will proceed by induction on A. First assume A = 1. Without any 
loss of generality (because of the scalar Casimir analogous to ( |3.5[ )) we may 
write 

lf> = £ C rs (J n °yV+r + J-Pl-I(j°n) + (J-) 2 Pk -2(Jl J~), 
r+s<k 

for some constants c rs . As J+ annihilates z n } acting with on z n ~ s , 
s = 0, . . . , k — 1, we find 

T « z n-s = (_i)« s i £ f ™Y Crs z « + i ower order terms, 

r=0 ^2/ 



leading to J2r=o(^Y c rs — f° r eacn s - It can be easily shown that the most 
general element in Span{( J°) r ( J+) s }^Zq satisfying this condition is 



Z r=0 

Finally, no monomial (J^) fc may be present in T^. Thus, 

T{ k) = q (n, l)p° h _ 1 (J+, J° n ) + ft(n, 1) ( P U(J° n ) + J-p k -2{Jl J')), 

completing the proof of (3.7) for A = 1. Now assume that (3.7) holds for 
every : V n — > V n -A- In particular, it is also true for an operator 
mapping V n into V n -A-i- As q a (n, A) (a = 0, . . . , A — 1) annihilates z n ~ s if 
< s < A — a, acting with T^+i on (s = 0, . . . , k) and reasoning as in 
the case A = 1, we deduce that a factor ( J° — n ~^+ Q ) must be present in front 
of each for a = 0, . . . , A. If we move these factors before their respective 
q a (n, A) and group together the monomials in £>fc_A-i(</°, J~) with no J - , 
we find that may be written in the form (3.7) with A replaced by A + l. 

ii). If we apply i) for A = k, we find = J2t=o°aQa(n, k), and acting 
on z n ~ s , s = 0, . . . , k, we conclude that c a = for all a. 

Hi). This is obvious, since derivatives of order higher than n annihilate 
V n . Q.E.D. 

The corresponding case : V n — > V ti +a is treated next: 
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Lemma 3.3 Let T^ k ' : V n — > V h+ a be a k-th order linear differential opera- 
tor. We then have: 

i) Ifn>k, then may be represented as 

Q = 

where each p k is a polynomial of degree not higher than k in the oper- 
ators J+, J° and J~. 

ii) If k > n, then T<® = T <9™ +1 + f, where T is a linear differential 
operator of order k — n — 1, and T is a linear differential operator of 
order at most n satisfying i). 



The proof of this Lemma is similar to that of Theorem [2.4| , and shall not 
be presented here. The Theorem analogous to |2.4| for a PVSP operator in 
now follows applying Theorem |2.4| and the Lemmas above to each entry 
of the operator: 

Theorem 3.4 Let n > m, and A = n — m. Let T^ k > be k-th order differential 
operator in 3^„. We then have: 

i) If m > k, then is a polynomial in the operators (3.1). More 
explicitly, if k > A > 1, 

T« = p fc (:n + Jp fc (T e ) + £ Q Q p a k _ A (T e ) + E QaP a k (T e ) , (3.8) 

o=0 o=0 

while if A = the same formula is valid with J replaced by J. Ifk < A, 
every p k _ A must be identically zero. 

ii) If n > k > m, then T^ k > = T d™ +1 + T, where T and T are matrix 
linear differential operators of the form 

/ a (fe-m-l) Q \ f~ a {m) ~ b (k)\ 

where the superscripts indicate the highest possible derivative in each 
entry, andT satisfies i). 
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Hi) If k > n, then = T <9™ +1 + T, where T is a 2 x 2 matrix linear 
differential operator of order k — n—1, and T : V m ^ n — > P m)n a linear 
PVSP operator of order at most n verifying i) or ii). 

The next issue to be addressed is to find out the number of parame- 
ters determining a generic fc-th order linear differential operator preserv- 
ing V m ,n-, that is, the dimension of 7^ n . In the scalar case, any k-th de- 
gree polynomial in J+, J°, J~ may be constructed from the monomials 
{(Jn Y(JnY~ r }r=oi s = 0, and so dim 3^ = (k + l) 2 if n > k, M 



Remarkably, in the matrix case we have dimT^ = A(k + l) 2 independently 
of m and n, provided m > k > n — m, as a consequence of the following 
Lemma: 



Lemma 3.5 The following set of monomials form a basis of the vector space 
of polynomials in the operators (3.1) of differential order at most k: 

{x{T ±y {TO y-r } s =Q ? {Q^)'}^ S = 0,...,k, 

where X = 1, J (or J, if A = 0), Qq, along with, if k > A: 

{QoO^r (T°r r };L , {Q a (T-) s }ti, s = 0, . . . , k - A. 

Proof. Linear independence of the monomials is straightforward from the 
definition of the operators. Completeness is a consequence of the following 
facts. In the first place, every J s is a linear combination of {1, J} (and 
analogously for J). Secondly, JQ a is proportional to Q a , and Q a Qp = 
(and the same for the Q's). Third, any product Q a Q/3 is a diagonal PVSP 
operator, and thus expressible through the T's and J (or J). Finally, the 
formulas (a > 1) 

QaT Qa—lT ~\~ ^ Qa > QaT Qa—lT ~\~ ~ Q a —\ , 

QaT = Q a _iT- + (- + l) Q a , Q a T + = Q a -iT° + ( - + l) Q a ^ , 

allow us to remove every T° and T~ (respectively T + ) from the monomials 
with Q a (respectively Q a ), a = 1, . . . , A. Q.E.D. 
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Corollary 3.6 Let n > m > k, and A = n — m. We then have: 
dirndl 



4(/t + l) 2 , k>A 
(Jfe + l)(3Jfe + A + 3), A>k. 



If m < k, dim '?^ n is no longer finite, as arbitrary differential operators are 
involved in this case. 



3.2 Case N > 2 

We now examine briefly some aspects of the case N > 2. Let n% < • • • < njv 
be non-negative integers, and let = rij — rij, where j > i. Consider the 
following set of N x N matrix differential operators, ||: 

T £ = diag(j: i ,...,^ iv ), e = +,0,-, 

Pi = diag(0, . . . ,0,1,0,. ..,0), 

Qa(i,j) = z a \ ij} i> j , a = 0, . . . , Ajj, 

Qa(i,j) = <la(nj, Ajj) Aij, j>«, Qf = 0, A^, (3.9) 

where (Ajj) pg = Si P 5j q . It can be readily verified that the operators in (3.9) 
preserve V ni ,...,n N - A complication arising when iV > 2 is to define a suitable 
composition law between the latter operators. In the approach of Brihaye 
et at, 0, this composition law is defined to be an anticommutator if both 
operators are off-diagonal and a commutator otherwise, but the algebra thus 
obtained is no longer a Lie superalgebra, since the anticommutator of two 
off-diagonal operators is not always a diagonal one. This reflects the fact 
that the Z 2 -grading we introduced for N = 2 (i.e. classifying the operators 
in diagonal and off-diagonal) does not define an associative superalgebra in 
D NxN when N > 2, for the usual product (composition) of two off-diagonal 
matrix differential operators is not necessarily diagonal. A possible general- 
ization of this Z 2 -grading, endowing D^xn with an associative superalgebra 
structure can be defined as follows. An operator T = aXij, where a e 2), is 
said to be even (respectively odd) if i + j is even (respectively odd). Hence, 
any diagonal operator is even. We can likewise use this grading and the gen- 
eralized Lie product ( |3.3| ) to construct a Lie superalgebra structure in T>n x n. 
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It is not clear, however, whether this construction is really useful, and so it 
will not be further discussed. 

As remarked by Brihaye et at, ||, Theorem |3.4| can be easily generalized 



to arbitrary N, the operators (3.9) playing the same role as those in (3.1) 
for N = 2. Moreover, it is not difficult to show that dim 3^ _ is still 
independent of the rij's if they are large enough and their differences are 
small enough. More precisely, if n\ > k, we have: 



dim 311.^ 



'N 2 (k + 1) 2 , k>A 1N 

N{N + 1} (k + l) 2 + + A 1N >k 
1 i<j 



where 8ij = Ajj if A^- > k, and 8ij = k + 1 if Ay < k. If k > m, arbitrary 
differential operators are involved and thus 7^ njv is infinite-dimensional. 
Although no attempt will be made here to give a formal definition of a 



QES algebra of matrix differential operators, it is clear that Definition [2.1| of a 
QES differential operator is too restrictive in the matrix case. Indeed, the re- 
sults of this section suggest that in the matrix case one should include at least 
Lie superalgebras of differential operators — not necessarily finite-dimensional 
nor spanned by first order operators — preserving a finite-dimensional mod- 
ule of functions among the class of matrix QES algebras. In any case, it 
is intuitively clear that PVSP operators are just a particular class of QES 
operators. 



4 Spin 1/2 Schrodinger Operators 

From now on we will deal only with 2x2 matrix second order differential 
operators (N = k = 2 in the notation of the previous sections). We start by 
formally defining the class of matrix Schrodinger operators: 

Definition 4.1 A Schrodinger-like operator is a second order differential 
operator of the form H = —d 2 + V(x), where V is an arbitrary 2x2 (complex) 
matrix. A Schrodinger operator (or Hamiltonianj is a Hermitian Schroding- 
er-like operator, i.e. the matrix V is of the form 

f Vl {x) v*(x)\ 

\v{x) v 2 {x)J ' l4 ' iJ 
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where Vi andv2 are real-valued functions andv is an arbitrary complex-valued 
function. 



The notion of equivalence we shall use for matrix differential operators 
is the same as in the scalar case (see Definition p.2| ), where now the gauge 



factor U(z) is an invertible complex 2x2 matrix. We will be interested in 
constructing one-dimensional Schrodinger operators H equivalent to a second 
order differential operator T in 7^ n , with A = n — m > 0. This equivalence 
can then be used to construct m + n + 2 eigenfunctions of H from the corre- 
sponding ones of T obtained by diagonalization of the (m+n+2) x (m+n+2) 
Hermitian matrix representing T in V mn . We will assume that m > 2, and 



thus T is a polynomial in the operators (3.1), according to Theorem CO 



Let T : V m ,n — ► V m ,n (n > m > 2) be a second-order PVSP operator. 



From Theorem !T3 we have 

— T — A 2 (z) dl + Ai(z) d z + A (z) , (4.2) 

where the A^s are 2x2 matrices with polynomial entries (in this section, 
capital letters will be reserved for matrices). Assume that T(z) is equivalent 
to a Schrodinger operator H(x) under a gauge transformation U(z) and a 
local change of variable by a real- valued function x = <p(z), i.e. 

-T(z) = -U- 1 (z)H(x)U(z) 

= d 2 x + 2Ad x -B, (4.3) 

with 

A(x) = U- 1 U x , B(x) = U~ 1 VU - A 2 - A X} U{x) = U{<p- 1 {x)) . 

Here and in what follows, a subscripted x denotes derivation with respect 
to x, while derivatives with respect to z will be denoted with a prime '. 
Expressing T(z) in the variable x, we obtain the operator T(x) given by 

- f(x) = \A 2V ' 2 d 2 x + (Anp' + A 2 v ") d x + A 1 u ^ , (4.4) 

L lip L (x) 

and comparing with (4.3) we conclude that A 2 must be a multiple of the 
identity. It then follows that only the first term in equation (3.8) of Theo- 
rem p.4| contributes to A 2 , and taking into account the explicit form of the 
T e 's (see (3.1)), we conclude that A 2 is a 4-th degree polynomial p 4 times 
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the identity matrix. (Unless otherwise stated, we will denote by p n {z) an 
arbitrary polynomial in z of degree at most n with complex coefficients.) We 
also deduce that (p(z) satisfies the equation p^' 2 = 1, or 

f z 1 

x = J ds } (4.5) 

/p 4 (s) 



and thus the coefficients of p 4 must be real. Identifying the corresponding 
remaining terms in (4.3) and (|4.4j) , we then get 



A(x) = ——(A 1 --j/A , B(x) = -A . (4.6) 



2^p-^ r 

Thus, we have shown: 



Theorem 4.2 Let T be a PVSP operator in J^ n , with n > m > 2. Then 
T is equivalent to a Schrodinger-like operator if and only if it is of the form 

-T = p 4 d 2 z +A 1 d z + A . (4.7) 

The operator T is equivalent to a Schrddinger operator —d 2 x + V(x) if and 
only if (|4.7|) holds, and in addition there is an invertible matrix U satisfying 
the differential equation 

U x = UA (4.8) 

and such that 

V = UWU~\ where W = B + A 2 + A x , (4.9) 
is Hermitian (with x, A, and B given by ( |4.5| ) and ( [4.6| ) ). 



The eigenfunctions of the Hamiltonian H are of the form i/?(x) = 
with ^(x) = \I/((y9~ 1 (x)), where ^(z) is an eigenfunction of T and ip must 
satisfy suitable boundary conditions to qualify as a physical wave function. 

We note that once an invertible solution U of equation ( |4.8| ) has been 
found, any other invertible solution is of the form UqII for some Uq in 
GL(2, C). In fact, multiplying U by such Uq is equivalent to performing 
a further constant scale transformation by Uq. In the scalar case this ad- 
ditional freedom is absent, for differential operators are unaffected by scale 
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transformations by constant functions. Note also that a scale transforma- 
tion by an arbitrary constant matrix Uq will not map every Hamiltonian into 
another Hamiltonian, unless Uq belongs to R + x U(2). 

A matrix differential operator will be called uncoupled if it is either upper 
or lower triangular. Since A2 in (|4.2| ) must be a multiple of the identity ma- 
trix, it follows that a PVSP operator of the form ( |4.7| ) will be automatically 
uncoupled whenever A > 1, as no Q can then be present in (3.8). Moreover, 
the following result shows that any Hamiltonian we may obtain from T when 
A > 1 will be essentially diagonal: 

Proposition 4.3 Every Hamiltonian H obtained from an uncoupled PVSP 
operator T of the form (|4.7| ) is diagonal, up to equivalence. 

Proof. If T is uncoupled, the integration of equation (|4.8| ) is straightforward. 
Multiplying U from the left by an appropriate Uq in M + x U(2), we construct 
a new gauge factor uncoupled in the same way as T. Using this new gauge 
factor, we obtain a Hamiltonian H given by 

H = UqHU, 1 = UqIIT^qUY 1 , 

which is both diagonal and equivalent to the initial one. Q.E.D. 

Consequently, we shall limit ourselves in what follows to the cases A = 
0,1. 

There are two main difficulties associated with the method just outlined 
for constructing QES spin 1/2 Hamiltonians. In the first place, one needs 
to invert the elliptic integral (|4.5|) in order to compute function of 

x, which is no easy task. Secondly, the differential equation ( |4.8| ) cannot 
in general be solved in closed form, thus preventing us from verifying the 
Hermiticity of V. The former complication can be overcome, as we shall see 
in the next section. The latter is more difficult to handle, although imposing 
further constraints on the initial PVSP operator will contribute to simplify 
the problem, as shown in Section |6|. 

We shall finish this section with a few remarks on the physical significance 
of matrix Schrodinger operators. First of all, one-dimensional 2x2 matrix 
Schrodinger operators can be obtained by separation of variables from the 
three-dimensional Pauli Hamiltonian describing a spin 1/2 charged particle 
in non-relativistic quantum mechanics. Consider, indeed, the Pauli Hamil- 
tonian 

#Pauii = (iV + e A) 2 + e - e cr • B , 
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where and A = (A 1 , A 2 , A 3 ) are respectively the scalar and vector potential 
of the external electromagnetic field, B = V x A is the magnetic field, the 
Pauli matrices cr = (a 1 , a 2 , a 3 ) are given by 

■■-(!;). H°~o')> 

e is the electric charge, and physical units have been chosen so that h = c = 
2m = 1. If, for example, the vector and scalar potentials depend only on the 
x coordinate (and we take, without loss of generality, A 1 = 0) then i?p au ii 
obviously commutes with the y and z components of the linear momentum. 
The eigenfunctions of i?p au ii can then be sought in the form 



e i{p y y+PzZ) if){x) 



where p y ,p z G M are the values of the y and z components of the linear 
momentum, and the two-component spinor if)(x) is an eigenfunction of the 
one-dimensional matrix Schrodinger operator with potential ( (4.1[ ) given by 

dA 2 

Vj ( x ) = e + (eA 2 - Py ) 2 + (eA 3 - p z ) 2 + (-l)'e — , J = 1,2, 
dA 3 

v(x) = ie — — . 

dx 

More surprisingly, one-dimensional 2x2 matrix operators are also directly 
related to Dirac's relativistic equation for a spin 1/2 charged particle in an 
external electromagnetic field. To see this, let us write the latter equation as 

(i# - e4 -m)*(a?) = 0, (4.10) 

where ^ = 7%^, the 7 M 's are 4x4 matrices satisfying 

the metric tensor {g^ v ) = diag(l, —1, —1, —1) is used to raise and lower in- 
dices, = d/dx^, x° = t, A = cf) and m is the particle's mass. Multiplying 
Dirac's equation by the operator i$ — eJf^ + m we easily arrive at the second- 
order equation 



(id - e A) 2 -m 2 - 6 -F^ 



tf = 0, (4.11) 
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where a^ v = | [7^, 7"] and = c^A, — d u A^ is the electromagnetic field 
strength tensor. Conversely, if $ is a solution of ([4.1 1|) then either — eJ^ + 
m)$ or 75$ = 27°7 1 7 2 7 3 $ (but not both simultaneously!) is a non-trivial 
solution of Dirac's equation (|4.10|) . In the chiral representation of the gamma 
matrices, ITTtl, we have 



a 



oj 











a 



. 







a 



(j, k,l = 1, 2, 3, and summation over I is understood), and therefore ( |4. 1 1|) de- 
couples into two independent equations for the upper and lower components 
V±{x) of namely (cf. @) 



[id - eA) 2 - m 2 + e a ■ (B =f % E) 



± 



0. 



If the electromagnetic four-potential is time-independent, and we look for 
solutions of Dirac's equation with well-defined energy E, i.e. we set \l/± = 
e~ lEt t{)±(x,y, z), we obtain the following equation for i/? ± : 



W + eA) 2 - (E 



eA ) 2 + m 2 - e a ■ (B =f i E)l if>± 







(4.12) 



which has the same structure as Pauli's non-relativistic equation. Just as was 
the case with Pauli's equation, separation of variables in ( |4.12[ ) often leads to 
the eigenvalue problem for a one-dimensional matrix Schrodinger operator. 
For instance, suppose that A = and that A has cylindrical symmetry, that 
is 

A = A{p) e z 

in cylindrical coordinates (p,(p,z). The left-hand side of ( |4.12|) then com- 
mutes with the z components of the linear momentum (—id z ) and the total 
angular momentum {—id v + |cr 3 ), which allows us to look for solutions of 
(gTj ) of the form 



iPzZ 



i2 2 (p)e^ +1 / 2 ^ 



(4.13) 



Here we have dropped the subscript ±, since if) + and satisfy the same 
equation, p z e R is the value of the z component of the linear momentum, 
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and j z G N + \ is the value of the z component of the total angular momen- 
tum. Substituting (|4.13| ) into (|4.12|) we arrive at the following equation for 



-cP - -d p + (p z - eAf + m 



+ 



E 1 
- 1/2) 2 











dA 2 
dp 



. R 2, 



Defining 



fi(x) =x l/2 Ri{x) 



1,2 



we finally obtain that the two-component spinor (fi(x) f2{x)^j is an eigen- 
function of the one-dimensional 2x2 matrix Schrodinger operator with po- 
tential 



V(x) 



Pz -eA(x)) 2 + (x) a 2 + ^ (j z -a 3 ), x G (0, oo) , 



with eigenvalue E 2 — m 2 and boundary condition /i(0) = /2(0) = 0. 



5 GL(2) Action and Canonical Forms 

In this section we will study how the GL(2) action on the projective line MP 1 
induces an automorphism in the superalgebras Sa generated by the opera- 
tors (3.1). This will allow us to reduce the polynomial Pa(z) to some simple 
canonical forms, facilitating the evaluation of the integral (|4.5| ). These ideas 
were first applied in the context of QES systems to analyze the normaliz- 
ability of the wave functions of scalar QES Hamiltonians, ||. We introduce 
some definitions and results in the scalar case, and then show how to extend 
these concepts to the matrix superalgebras Sa- 

The action of GL(2) = GL(2, R) on MP 1 via linear fractional or Mobius 
transformations, 




(5.1) 



induces an action on V n , mapping a polynomial p(w) to the polynomial p(z) 
given by 

p( z ) = ( 7Z + S)"p(?LtP y ). (5.2) 
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This defines an irreducible multiplier representation of GL(2) in ?„, [Id" 



which will be denoted by p n $. Since the infinitesimal generators of this mul- 
tiplier representation coincide with the generators of Q n , cf. fl2.4|) , it follows 
that the representation p n induces an automorphism of the Lie algebra g n 
spanned by the J^'s in Q2.5|) , ||. Performing the explicit scale transformation 
and change of variable, 



we obtain: 



1 

c\ 



yz + 5 . 

fa 2 2ap 
\ 7 2 2 7 5 



(5.3) 




Therefore, the J^'s transform according to the representation P2-1 = P2,o ® 
det -1 of GL(2), where det -1 is the reciprocal of the representation det : C 1— > 
|C|. It is convenient at this stage to introduce a larger class of representations 
ofGL(2): 

Definition 5.1 Let n > 0, i be integers. The (irreducible) multiplier repre- 
sentation p n j of GL(2) on V n is defined by 



p{w) h-> p(z) = (a5 - f3y) % (jz + S) n p(^ 



az + (3 
■yz + 5 



We note the isomorphism between p nyi and p n ® det*. As shown in ||, a 
second-degree polynomial in the J^'s (in fact, any operator in 7^ if n > 2, or 



any QES operator on the line modulo equivalence, according to Theorems |2T3 
and p.4|) may be written as 



V 



/■?. 



77 



— p'J <9 Z + r - -g' + 



n(n 



12 



(5.4) 



where p, q, and r are polynomials in z of degrees 4, 2, and respectively. The 
transformation of P2{J^) under the action Q5.3|) is easily described in terms 
of the triple (p, q, r), [|]: 

Lemma 5.2 Lei P2 be a second-degree polynomial in the operators J^(w), 
determined by the triple (p(w) , q(w) , r) . Then, the transformed polynomial 
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P2 under the GL(2) action ( |5.3| ) is determined by the triple (p(z) , q(z) , r) 
given by 

(>yz + 5) 4 faz + l3\ (jz + 5) 2 faz + (3 



Piz) = ^c\^ p {^zT5)> ~ q{z) = ^cT q K^5)> r = r - 

Therefore, a second-degree polynomial P2 in the J^'s transforms accord- 
ing to the direct sum representation P4,_2 © P2-1 © Po,o under the GL(2) ac- 
tion (5^3). One can choose a particularly simple representative of the GL(2) 
orbit generated by P2 by placing the polynomial p (assumed to be real) in its 
associated triple (p, q, r) in canonical form, 



Theorem 5.3 Every non-zero quartic real polynomial p(z) transforming un- 
der the representation p±-i of GL(2) is equivalent to one of the following 
canonical forms: 

1) v(z A + tz 2 + 1) , r ^ ±2, 5)z/(^ 2 -l), 

2) u(z 4 + tz 2 - 1) , 6) vz 2 , 

3) u(z 2 + l) 2 , 7)z, 

A)u(z 2 + 1), 8)1, (5.5) 

where v 7^ and r are real numbers. 

We now generalize these results to the matrix case. The induced action 
of GL(2) on V n -A,n analogous to ( |5.2| ) is: 



yP \w) VP\z) v 1 \p 2 (w(z)) / 



where 

U(z) = diag(( 7 z + 5) n - A ,( 1Z + 5) n ), 

and w(z) is given by ( |5.1| ). This representation of GL(2) in V n -.A,n is obvi- 
ously isomorphic to p n -A,o©Pn,o- The following Lemma describes the induced 
action on 5a: 

Lemma 5.4 The action o/GL(2) on Sa given by 

X(w) h- X(z) = £>(*) U-\z) , X G s A , (5.6) 
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defines a Lie superalgebra automorphism. The generators of s&, cf. (3.1) 7 
transform according to the following irreducible representations : 



{T 6 } -> 02,-1 , U}^Po,o, 
{Qa} —> PA,0 , {Qa} —> PA, -A 



A straightforward generalization of equation ( |5.4j ) and of Lemma p.2\ to 



a second-degree polynomial in the T e 's shows that the (real) polynomial p^ 



in (|4.7| ) transforms according to the representation Pa- 2 under the GL(2) 
action ( |5.6| ). We will henceforth assume that p± is one of the canonical forms 
given in Theorem |5.3| . The integral Q4.5|) and the inverse z = y? _1 (x) can 



then be easily computed for each of these canonical forms, . 

Before finishing this section, let us point out that equations Q4.8D and ( |4.9| ) 
adopt a simpler form in the variable z, because in that case only rational 
functions appear. Explicitly, the equation for the gauge factor reads: 

U' = UA, with A(z) = ^ = 7 ^-(A 1 -U), (5.7) 



while U and W in equation (|4.9|) are substituted by U and W, where 

W = -A + p 4 A' + p 4 A 2 + ^p^A. (5.8) 



We now use Theorem |3.4|, Lemma 1313 and the explicit form of the oper 



ators (3.1) to compute A and A for the most general operator in !P^ A(| 



of the form (|4.7|) , in the cases A = 0,1. We denote by p" the polynomial 
J2?=o a i z \ where on are arbitrary complex numbers. If the polynomial P4 is 
not one of the first three canonical forms, we obtain: 



Case A = 0: 



Pi A 



P2 V\ 



K P 2 Pi, 
Case A = 1: 



Pa A 



( P\ P b i 

yP3 Pi, 



An 



A 



do — 2na 2 z b — 2nb 2 z 
c — 2nc 2 z do — 2nd 2 z, 



do — 2(n — l)a 2 z 
c + c.\Z — 2(n — l)c 3 z 2 



(5.9) 



-2nh N 
do — 2nd 2 z 



(5.10) 
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where a ,b ,c ,Ci, and d are arbitrary complex numbers. If p 4 is one of the 
first three canonical forms, the following extra terms are present in p^A and 

p A A -> -diag (jji - A)uz 3 ,nuz 3 ^j , (5.11) 
A -> + diag ((n - A) (n — A - 1)zaz 2 , n(n - l)z/z 2 ) . (5.12) 



6 The Gauge Factor 

In this section we will deal with the differential equation (|5.7|) for the gauge 
factor U(z). As remarked in Section 0, this equation cannot be solved in 



closed form for every A. Equation (|5.7|) admits the formal iterative solution 

/ J^L r z r s i rs„-i A \ 

U(z) = U (l + YlJ ds i J ds 2 ---J ds n A(s n )...A(s 2 )A( Sl )j, 

but this is of little practical use in checking the Hermiticity of V. We will 
impose a certain condition on A so that Q5.7T) can be explicitly solved, thus 
reducing the number of parameters defining A, but still leaving us with room 
to construct relevant examples of QES Schrodinger operators. 
Denoting 

u=h u A , A=( a b )i , 

we may write equation ( |5.7|) as 




(6.1) 

with the same system for (7x3,1*4)*. Unfortunately, the associated scalar sec- 
ond order differential equations for u\ or w 2 are not any of the standard equa- 
tions of Mathematical Physics. We also note that the quotient q = U1/U2 
satisfies the following Riccati equation: 

q' = -bq 2 + (a-d)q + c. (6.2) 



In fact, solving (|6.1|) is equivalent to solving ( |6.2|) . Moreover, if we try to 
uncouple A by performing the linear transformation 
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we find that the transformed of A will be lower triangular whenever q satisfies 



qs' — sq' + cs 2 — bq 2 + qs(a — d) = , (6.3) 

with a similar condition to change A into upper triangular form. If s is non- 
zero we can take s = 1 without any loss of generality, and equation ( |6.3|) 
reduces to ( |6.2| ). Unfortunately, none of these equations can be solved with- 
out further assumptions. 

However, if we restrict ourselves to matrices A satisfying the equation 

[A(z), T A(s)ds] = Q, (6.4) 
for some z G R, we can readily integrate the gauge equation (|5.7|). Recall 



that this condition on A was indeed verified by the QES Schrodinger operator 
found by Shifman and Turbiner, [T7|]. If (|6.4j) is satisfied, we shall say that 



we are in the commuting case. In this case, the general solution of the gauge 
equation is given by: 

U(z) = U exp f A(s) ds , (6.5) 
where Uq is in GL(2, C). It is worth mentioning in passing that we can use 

£>0, r/>0, 0g[O,2tt) (6.6) 



' 1 X 



to parametrize the orbits of R + x U(2) acting on GL(2, C) by left multi- 
plication. Since matrices in the same orbit lead to equivalent Schrodinger 
operators, we can take Uq in the form ( |6.6| ). 

We shall be mainly concerned with the commuting case. We will make 
use of the following elementary Lemma to describe the most general form of 
A in the commuting case: 



Lemma 6.1 Let M(z) be a 2x2 matrix satisfying the equation ( |6.4|) . Then 
M is of the form: 

M = f{z)M +g(z), (6.7) 
where f and g are scalar functions, and M is a 2 x 2 constant matrix. 
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With this Lemma in mind, looking at the expressions for A (A = 0, 1) 
that we obtained in Section [5[ cf. (5.9)— (5.11), we find that the most general 
A satisfying (|6.4j ) is of the form 

Pi A = p 2 (z) + A(z), 

where, as usual, pi denotes a second-degree polynomial in z with complex 
coefficients, and the matrix A is given in Table 1. 

P4 in canonical forms 1-3 p^ in canonical forms 4-8 



a 0\ 

7 S) 

'a (A [0 (A 

7 fj'"^ Oj 

Table 1: Matrix A(z) for the canonical forms 1-8 (at, /3, 7, S G C). 

Note that if A = 1 and p^ is one of the first three canonical forms, every 
Hamiltonian we can possibly obtain will be diagonal, modulo equivalence. 

Finally, let us remark that in the non-commuting case (that is, when 
[A, J z A] 7^ 0), we may still be able to integrate Q5.71) explicitly by imposing 
other constraints on A, as e.g. assuming it is uncoupled. Alternatively, if p^ is 
not any of the first three canonical forms, and we assume that the columns (or 
rows) of A are proportional to each other (the ratio of the respective entries 
being a constant), we can also reduce (|5.7|) to quadratures. Unfortunately, 
we have not been able to find any interesting examples of QES Hamiltonians 
in the non-commuting case. 



An ( a P\ 

A = Do _ — nvz'- 

V7 



P2 



A = 1 



(p 2 - isz 3 



n-1 

i 7 n. 



a {3\ 



7 Examples 

In this final section we exhibit some new examples of spin 1/2 Schrodinger 
operators equivalent to a PVSP operator of the form (JO). In the previous 
section we have seen how, by restricting ourselves to the commuting case, we 
were able to integrate equation ( |5.7| ) explicitly. This is not, however, the end 
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of the problem, for we must still check that the matrix 



V = UWU' 



with W given by (|5.8|), is self-adjoint. Again, if we start with the most 
general PVSP operator of the form ( |4.7| ), the algebraic constraints imposed 
by the condition V — are too complicated to be solved in full generality, 
even if we limit ourselves to the commuting case. This situation is completely 
analogous to what we find when trying to construct scalar QES Schrodinger 
operators in more than one spatial dimension, [[HJ, |7]]. We have thus no 
choice but to look for particular examples. We now present some relevant 
examples of QES spin 1/2 Hamiltonians for the commuting case. In the first 
two examples we take A = 0, while in the remaining ones we assume that 
A = 1. Finally, let us point out that many otherwise interesting examples 
are often reduced to trivial ones after we impose the square integrability con- 
dition on the eigenfunctions. 

Example 1. Consider the six-parameter PVSP operator T given by 

-T = T-T° + (a 2 + d 2 ) T + + (a 2 - d 2 ) JT + + 2a x T° + (2a Q + ~(n - 1)) T~ 

+ (Qo + Qo)(2 6 2 T + + 6 ), 
where all the parameters are real numbers coinciding with those appearing 



in equation (5.9), and b 2 ^ 0. Note that we are in case 7 of Theorem 5^ 
(p4 = z), and so z = x 2 /4. The gauge factor U(z) is given by 

rrf \ — f< \ ( A cos h^ + (ct2 — d 2 ) sinhw 2b 2 sinhw \ 

\ 2b 2 sinhu A coshu + (d 2 — a 2 ) sinhuy ' 

with 

f(z) = z ao exp (aiz+j(a 2 + d 2 )z 2 ) , u = jz 2 , A 2 = (a 2 -d 2 ) 2 +4b 



2 

2 ■ 

where we have taken U = A in equation ( |6.5| ). Using ( |7. 1| ) , we obtain a 
potential V(x) with entries given by (see (|4.1| )): 

v\ = a^ 2 x~ 2 + «o + Oi 2 x 2 + a^x 4 + a^x 6 , 
v 2 = 5- 2 x~ 2 + 5 + 5 2 x 2 + 5 4 x 4 + 5 6 x G , 
v = fa + p 2 x 2 + p 4 x 4 + (3 6 x 6 , 
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where 



a 2 bo 0^02/ ,3^ 

a\a 2 a\ + b\ 



a_ 2 = 2a (2a - 1) , a = — — , a 2 = — + — (a + n + 



« ? o 62/ , , 3x _ b 2 (a 2 + d 2 ) 
Po = -b , ^ 2 = -(a + n + i j, ^4 = -g- , ^ 6 = — , 

while each S n is given by the corresponding a n replacing a 2 by d 2 . We have 
ignored a constant multiple of the identity in V, which is equivalent to fixing 
a new origin in the energy scale (this will be done for the subsequent examples 
without further notice). 

We note that the potential V is singular at x = unless a is either or 
1/2. Let us introduce the parameter A = 2ao — 1, in terms of which we have 
a_ 2 = A(A + 1). If A is a non-negative integer /, we may regard 

{-d% + V(x)-E)ij)(x) = 0, 0<x<oo, (7.2) 

as the radial equation obtained after separating variables in the three-dimen- 
sional Schrodinger equation with a spherically symmetric Hamiltonian given 
by 

H = -A + U(r) , with U(x) = V(x) - fcil , 

x 2 

where A denotes the usual flat Laplacian. Given a non-negative integer I and 
a spherical harmonic Yi m (9, 0), if ip is an eigenfunction for the equation ( |7.2| ) 
satisfying 

lim ij>(x) = 0, (7.3) 

x— »0+ 



then 



*(r,e,<t>) = ^-Y, 



lm \ 



will be an eigenfunction for H with angular momentum I. If A is not a 
non-negative integer, we shall consider (|7.2p as the radial equation for the 



singular potential U(r) = V(r) at zero angular momentum. The potential 
U(r) is physically meaningful, in the sense that the Hamiltonian H admits 
self-adjoint extensions and its spectrum is bounded from below, whenever 
A 7^ —1/2, ||, 0. The boundary condition Q7.3Q must be satisfied in the 
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singular case for all values of A. This boundary condition is verified if and 
only if a > 0. Finally, the additional conditions 

a 2 d 2 > &2 5 a 2 < , 

ensure that xj) lies in L 2 (I) © L 2 (I), where I = [0, oo) in the singular case 
or at zero angular momentum, or / = R in the non-singular one-dimensional 
case. 



Example 2. As our second example, we take 

-T = (T~) 2 + (ai + di) T° + ( ai - d x ) JT° + 2a T~ 

+ (61Q0 + blQ )(2T° + n -fi) + ~{n- fx)( ai - d x )J, 

where the coefficients are real, excepting b\ ^ which is complex. In this 
example p± — 1 and thus z = x (case 8 of Theorem jOf ). We choose the 
following gauge factor: 



A coshu + (ai — c?i) sinhw 2 6^ sinhu N 

2 6x sinhw A coshw + (d\ — a x ) sinhu 



U{z) = f{z) 
with 

f(z)=exp(a z+^(a 1 + d 1 )z 2 ^, u = ^z 2 , A 2 = (ai - c/i) 2 + 4 |&i| 2 . 

We obtain the potential with entries given by 

vi = ai([i + 1) + 2a a 1 x + (a 2 + |&i| 2 ) x 2 , 
^2 = di(fi + 1) + 2a d 1 x + (d\ + |&i| 2 ) x 2 , 
v = bx (u + 1 + 2ao x + (ai + di) x 2 ) . 

The eigenfunctions ip(x) are square-integrable provided 

a\ di > \bi\ 2 , ai < . 

Since the parameter n does not appear in the gauge factor or in the poten- 
tial, the Hamiltonian preserves U ■ V n , n for arbitrary n. (In the literature, 



an operator with this property is usually referred to as exactly solvable, |19 
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3D], [1]].) Note that if fi = 2n we could also obtain T from the PVSP op- 
erators (3.1) with A = 1, but in this case the Hamiltonian would preserve 
U ■ V m -i,m om y f° r rn = n. 

Example 3. From now on we take A = 1. Let T be the four-parameter 
PVSP operator given by 

-T = (T ) 2 + 2a 2 T+ + 2(n + 1)T° - 2 JT° + 26j Q + 26 Qi - 26j Q T° 
- 2b Q T- - (Aa 2 b + (3n + 1)6 X ) Qo - 4a 2 &iQi - (2d + n + -) J, 



with all the parameters real. Since = z 2 we are in case 6 of Theorem |573 . 
Solving equation ( |4.5| ) for z, we obtain z = e x . The gauge factor reads: 

tti \ r aoz ( cosu sinu\ x b 

U(z) = Jze 2 . , where u = \-bx\ogz. 

l-sinu cosu J z 

The potential is given by 

v . = -b 2 Q e - 2x - 2b h e~ x + (2n + l)a 2 e x + a\ e 2x 

+ (a(x) cos2w - (3(x)sm2uj , j = 1,2 

u = a(x) sin 2-u + COS2-U, 

where u = b\ x — b e~ x , and 

"0*0 = -y + 026^, /3(x) = (2n + 1) 6i + 2a 2 (6 + &i e x ) . 

It may be easily verified that the expected value of the potential is bounded 
from below, i.e. 

(ip,Vij>) > c\\*l>\\ 2 , with ip e L 2 (M)©L 2 (R), (7.4) 

for some c 6 1, if and only if fe = 0. (Note, however, that even in this 
case the amplitude of the oscillations of v(x) tends to infinity as x —>■ +oo.) 
Finally, the condition a 2 < is necessary and sufficient to ensure the square 
integrability of the eigenfunctions ij)(x). 
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Example 4. As our last example, we consider: 

-T = T-T° + 2 ai T° + (2a + n - ~) T~ - JT~ + 2h Q - 2b t Q T° 
- 61 (4a + 3n + 1) Q - Aa 1 b 1 Q 1 + 2(2a - a x ) J , 

where all the coefficients are real, and b\ 7^ 0. Since p 4 = z (case 7), we have 
z = x 2 /4. The gauge factor is chosen as follows: 

u{z) = zao eai J coshz sinM\ 
y— sinoi-z costly 

The entries of the potential V(x) are given by 

2a (2a -l) 1, 2 2 , 2 , & i x2 / \ • ^ 
= 2 !~ 7 ( a i ~ "1) x + ( °o cos —5 ot(x) sin — — - 

f = a sin — — + cos — — , 

with j = 1, 2, and a(x) is defined by 

a(x) = — (4a + An + 1 + aix 2 ) . 

We first note that the potential is singular at the origin unless = 0, 1/2. 
This situation is completely analogous to that of Example 1, so it will not 
be further discussed. The expected value of the potential is bounded from 
below, that is, equation Q7.4Q holds, if and only if 



a 1 



>1 + V2. 



Finally, the conditions 

a > , ai < 



guarantee the square integrability of the eigenfunctions if)(x) 
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